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Abstract 

The transition energy from galactic to extragalactic cosmic ray sources is still un- 
certain, but it should be associated either with the region of the spectrum known 
as the second knee or with the ankle. The baseline design of the Pierre Auger Ob- 
servatory was optimized for the highest energies. The surface array is fully efficient 
above 3 x 10^*^ eV and, even if the hybrid mode can extend this range below 10^^ 
eV, the second knee and a considerable portion of the wide ankle structure are left 
outside its operating range. Therefore, in order to encompass these spectral fea- 
tures and gain further insight into the cosmic ray composition variation along the 
transition region, enhancements to the surface and fluorescence components of the 
baseline design are being implemented that will lower the full efficiency regime of 
the Observatory down to ~ 10^^ eV. The surface enhancements consist of a graded 
infilled area of standard Auger water Cherenkov detectors deployed in two triangu- 
lar grids of 433 m and 750 m of spacing. Each surface station inside this area will 
have an associated muon counter detector. The fluorescence enhancement, on the 
other hand, consists of three additional fluorescence telescopes with higher elevation 
angle (30° — 58°) than the ones in operation at present. The aim of this paper is 
threefold. We study the effect of the segmentation of the muon counters and find 
an analytical expression to correct for the under counting due to muon pile-up. We 
also present a detailed method to reconstruct the muon lateral distribution function 
for the 750 m spacing array. Finally, we study the mass discrimination potential of 
a new parameter, the number of muons at 600 m from the shower axis, obtained by 
fitting the muon data with the above mentioned reconstruction method. 
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1 Introduction 

The cosmic rays energy spectrum extends for about eleven orders of magni- 
tude, starting at energies below 1 GeV up to energies of more than 10^° eV. It 
presents three main features: the knee, the second knee and the ankle. There 
is evidence of a fourth feature, the so-called GZK suppression [1,2], which is 
originated by the interaction of high energy protons with the photons of the 
cosmic microwave background. In the case of heavier nuclei, a similar effect is 
expected due to the fragmentation of the nuclei in their interaction with the 
photons of the microwave and infrared backgrounds [3] . 

The knee has been observed by several experiments [4,5,6] at around 3 — 5 x 
10^^ eV. At this energy the spectral index changes from —2.7 to —3.1. The 
KASCADE data shows that the composition at the knee presents a transition 
from light to heavy primaries in such a way that, at energies above 10^^ eV, the 
composition is dominated by heavy nuclei. These particles are originated in 
our Galaxy and what is being detected is, very likely, the end of the efficiency 
of supernova remnant shock waves as accelerators. 

The second knee has been observed at around 4 x 10^'' eV by Akeno [7], Fly's 
Eye stereo [8], Yakutsk [9] and HiRes [10]. The physics of this feature is still 
unknown, it might be due to the end of the efficiency of supernova remnant 
shock waves as accelerators or a change in the diffusion regime in our Galaxy 
[11,12]. 

The ankle is a broader feature that has been observed by Fly's Eye [8] and 
Haverah Park [13] , centered at approximately the same energy, ~ 3 x 10^^ eV. 
These results have been confirmed by Yakutsk [9], HiRes [10] and Auger in 
Hybrid mode [1]. AGASA also observed the ankle but at higher energy, around 
10^^ eV [14]. As in the case of the second knee, the origin of the ankle is still 
unknown and its physical interpretation is intimately related to the nature 
of the former. The ankle could be the transition between the Galactic and 
extragalactic components [15] or the result of pair creation by cxtragalactic 
protons in the interaction with the cosmic microwave background [16]. 

The precise determination of the mean chemical composition of the cosmic 
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rays in the energy range above ~ 10^ eV will allow us to understand the 
origin of the second knee and the ankle and to know the energy and the speed 
at which the transition between the Galactic and extragalactic components is 
given [17]. In particular, it will permit to decide among the three main models: 
(i) the mixed composition scenario [15], in which the composition injected by 
the extragalactic sources is assumed to be similar to the one of the Galactic 
sources and in which the transition takes place in the ankle region, (ii) the dip 
model [16] , in which the ankle is originated by the interaction of extragalactic 
protons with the cosmic microwave background and the transition is given at 
the second knee and (iii) a two-component transition from Galactic iron nuclei 
to extragalactic protons, around the ankle energy [18]. 

The Pierre Auger Observatory consists of two Observatories situated one in 
each hemisphere. The Southern Observatory, located in Pampa Amarilla close 
to the city of Malargiie, Province of Mendoza, Argentina, currently consists of 
nearly 1600 Cherenkov detectors placed in a 1500 m triangular grid covering an 
area of 3000 km^ plus four fluorescence telescope buildings, with six telescopes 
each, situated in the periphery of the surface array and overlooking it. The 
construction started in 2000 and is going to be completed early in 2008. A 
complementary Northern Observatory will be sited in Colorado, United States 
of America. 

The Southern Observatory, in its original design, is able to measure cosmic 
rays of energies above 3 x 10^^ eV for the surface array and < 10^^ eV in 
hybrid mode. Two enhancements, AMIGA (Auger Muons and Infill for the 
Ground Array) [19] and HEAT (High Elevation Auger Telescopes) [20], will 
extend the energy range down to lO^'^ eV, encompassing the second knee and 
ankle region where the Galactic-extragalactic transition takes place. 

AMIGA will consist of 85 pairs of Cherenkov detectors and 30 m^ muon coun- 
ters buried ~ 2.5 m underground, placed in a graded infill of 433 m and 750 
m triangular grids. The AMIGA infill area is bound by two hexagons cover- 
ing areas of 5.9 km^ and 23.5 km^ corresponding to the 433 m and 750 m 
arrays, respectively. The energy thresholds of the 433 m and 750 m arrays are 
~ 10^^ eV and ~ 10^^-^ eV, respectively [21]. On the other hand, HEAT will 
be formed by three additional telescopes of 30° to 58° elevation angle located 
next to the fluorescence telescopes building at Coihueco. They will be used in 
combination with the existing 3° to 30° elevation angle telescopes at Coihueco 
as well as in hybrid mode with the AMIGA infills. 

These enhancements will also allow detailed composition studies based on 

the combined measurement of the atmospheric depth of maximum shower 
development, Xmax, and the shower muon content. These two parameters are 
very sensitive to primary mass composition. Other mass sensitive parameters, 
like the slope of the lateral distribution function, rise- time of the signals in 
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the surface detectors, curvature radius, etc. strongly depend on them [22]. 



In this paper we will concentrate on the AMIGA nmon detectors [19]. These 
counters will consist of highly segmented scintillators with optical fibers end- 
ing on 64-pixel multi-anode photomultiplier tubes (PMT). The scintillator 
strips will be equal to those used for the MINOS experiment [23] . The current 
basehne design calls for 400 cm long x 4.1 cm wide x 1.0 cm high strips of 
extruded polystyrene doped with fluors, POP (1%) and POPOP (0.03%), and 
co-extruded with Ti02 refiecting coating. They are covered with refiective Al 
foil. To extract the scintillation light, a wavelength shifting fiber is glued into 
a grove which is machined along one face of the scintillator strip. A 10 m^ 
module will consist of 64 scintillator strips with the fibers ending on an opti- 
cal connector matched to a 64-pixcI multi-anode Hamamatsu H7546B PMT 
of 2 mm X 2 mm pixel size, protected by a PVC casing. Each muon counter 
will consist of three of these 64-channel modules, totalling 192 independent 
channels covering an effective area of 30 m^ (actually during the engineering 
array phase one of these 10 m^ modules from each counter will be split into 
two 5 m^ modules for further analyses close to the shower core). These muon 
counters will be buried alongside a water Cherenkov tank. Each of the 192 
channels of the electronics will count pulses above a given threshold, with an 
overall counter time resolution of 20 ns. 



We also present a detailed method for the reconstruction of the Muon Lateral 
Distribution Function (MLDF) from data obtained by the muon counters of 
the 750 m-array. An associated problem is the pile-up effect due to the finite 
segmentation of the muon counters. We analyse this problem and propose a 
correction that considerably improves the reconstruction of the MLDF. The 
number of muons at 600 m from the shower axis, A^^^'^(600), is extracted 
from MLDF fits using our reconstruction method. Subsequently, the design 
parameters of the muon counters (segmentation and area) arc validated by 
studying the impact of these parameters on the total A'^^'^(600) uncertainty. 
Finally, we study the mass discrimination power of A^^^^(600) as compared 
to other parameters normally used in composition analyses: the maximum 
development of the longitudinal profile, X^^ax, the curvature radius of the 
shower front, R, rise time of the signal at the water Cherenkov detectors, ti/2, 
and slope of the lateral distribution function of the total signal deposited in 
the water Cherenkov detectors, f3. Second order effects, like multiple triggering 
due to electrons scattered by a single muon, will be dealt with in a subsequent 
work. 
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2 Muon Counter Segmentation 



The AMIGA counter electronics just counts pulses above a given threshold, 
without a detailed study of signal structure or peak intensity. This method is 
very sturdy since it docs not rely on deconvoluting the number of muons from 
an integrated signal. It does not depend on the PMT gain or gain fluctuations 
nor on the muon hitting position on the scintillator strip and the corresponding 
light attenuation along the fiber track. Neither does it require thick scintil- 
lators to control Poisson fiuctuations in the number of photo electrons per 
impinging muon. But this one-bit electronics design relies on a fine counter 
segmentation to prevent undercounting due to simultaneous muon arrivals. 



2.1 Characteristic distance to the nearest station 

The optimum segmentation of the muon counters depends on the number of 
muons arriving in the time given by the system time resolution. Since the 
MLDF decreases very rapidly with the distance to the shower axis we need 
to determine the average position of the closest station, where pile-up is more 
significant. 

The position of the first station depends on the geometry of the array and 
the angular distribution of cosmic rays. We performed a simple Monte Carlo 
calculation for the 750 m-array. We uniformly distributed impact points in a 
triangular grid with arrival directions following a sin(^) x cos(^) distribution 
with zenith angle 9 < 60° and a uniform distribution for the azimuth angle <p. 
For each event we obtained the distance (along the shower plane) between the 
shower core and the closest station. Figure 1 shows the obtained distribution, 
with a mean value at ~ 230 m. Therefore, we will base our subsequent studies 
on the number of muons found at a characteristic distance of 200 m. 



2.2 Underground muons at the nearest station 

To prevent contamination due to charged electromagnetic particles of the 
shower, the muon counters will be buried underground. Muons lose energy 
mainly by ionization when they propagate through the soil. Assuming that 
the energy loss is proportional to the track length, the energy of a muon that 
travelled a distance x is given by, 

E^{x) = E^Q - apx, (1) 
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where E^o is the initial energy of the muon, we assume a soil density oi p — 1.8 
g cm~^ and a = 1.808 MeV cm^ g^^ is the fractional energy loss per depth of 
standard rock {{Z) = 11 and (A/Z) = 1/2) [24]. 

Showers initiated by iron nuclei produce more muons than those from lighter 
nuclei of the same energy. For that reason the simulations were performed with 
iron showers to take into account the most unfavorable case. We assumed 20 
ns of sampling time for the system time resolution. 

Figure 2 shows the time distribution of muons in bins of 20 ns at 200 m from 
the shower core, in an area of 30 m^ x cos 9 (muon counter area projected 
on the shower plane) and at 2.5 m depth for simulated iron showers of ^ = 
30°, primary energy E = 10^^ eV and using QGSJET-II [29,30] as the high 
energy hadronic interaction model. We used Aires version 2.8.2 [31] to simulate 
air showers and, to obtain the muon time distribution, we propagated them 
through the soil assuming that they move at the speed of light and the energy 
loss is given by Eq. (1). Note that only muons with E^q > O.SGeV/ cos^j {9i 
zenith angle of the individual muons) reach the detector. 

From figure 2 we can see that the maximum number of muons at ii^ = 10^^ 
eV and 9 — 30° in the first 20 ns-bin is about 90. 
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Fig. 2. Time distribution of muons at 2.5 m underground and at 200 m from the 

shower axis in an area of 30 x cos (30°) m2. The showers were initiated by iron nuclei 
oi E = lO^s eV and 6* = 30°. We used Aires with QGSJET-II to simulate the 50 
showers. 



2.3 Muon pile-up correction 



The probability distribution corresponding to a given number of muons that 
hit a segmented detector in a time interval is multinomial, 



P(ni,...,njv,,J 



N„ 



nil . ..um 



(2) 



where Ngeg is the number of segments, A^^ is the total number of muons and 



rii is the number of muons that hit the i— th segment, i. e. X^^i" rii — Nf 



As the AMIGA muon counter electronics will be designed to only count events 
above a certain threshold, without analysing the corresponding signal, when- 
ever two or more muons hit the same scintillator strip in the same time bin, 
these multiple muons will be counted as one. Therefore, the total number of 
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muons counted is given by, 



Ns. 



(3) 



i=l 



where Q{x) = if a; = and Q{x) = 1 if a; > 1. 

By using Eq. (2) we can calculate the mean value and the variance of , 



<n;^>{n,)=k 



seg 



N. 



seg, 



Var 



n: 



c 



seg , 



seg , 



(4) 



1 + {Nseg - 1) 1 



TV — 1 

^^seg 



(5) 



Figure 3 shows the number of counted muons as a function of the incident 
muons for a muon counter with 192 segments. 

To correct the effect of segmentation we can invert Eq. (4) , 



Corr 



In 1 - 



Nsea 



(6) 



The number of muons inferred, N^°'^'^, obtained from Eq. (6) has an error 
which can be estimated by solving for a given the following equations. 



(7) 
(8) 



where iV^*" and N^"-^, the unknowns, are the upper and lower limits of the 
error interval, respectively. The error associated with N^'°^^ is asymmetric, this 
is due to the shape of as a function of N^. If we define (j+ = N^"'^ — N^°'^'^ 

and (7_ = N'^'^'^ — A^,, we can see that a+ > (T_ and that the difference 

M ' ' 

increases with A^,, . 

To assess the importance of the error introduced with the use of the correction 
formula, it has to be compared to the Poissonian fluctuations inherent to the 
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Fig. 3. Mean value and one sigma region of the number of muons counted as a 
function of the incidents for 192 segments (three PMTs with 64 pixels each). In the 
ideal case of infinite segmentation = N^^. 



process of counting the muons that hit each scintillator strip in a given time 
bin. Therefore, the total error in the determination of the number of muons 
is AA^"^ = y'cr^ + (Jp^-gg, where crpoiss = is the error corresponding to 

the Poissonian fluctuations. 



Figure 4 shows the ratio between the total error and the Poissonian one as a 
function of the number of muons obtained after correction for 192 segments. 
From this figure we see that for 90 incident muons the total error is greater 
than the Poissonian by less than ~ 14%, and as such there is no need to 
further segment the detector. Note that the individual uncertainty in each 
muon counter will not directly translate to the iV^(600) but that it will be 
further reduced when a number of counters are used to fit the MLDF. 



It is worth emphasizing that AMIGA envisages an experimental verification 
of the segmentation, since in its unitary cell with 7 counters, the modules will 
have double segmentation. 
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Fig. 4. Ratio between the total error and the Poissonian error as a function of the 
corrected number of muons, for a counter with 192 segments. The two curves show 
exphcitly the asymmetry of the errors -see the text for further details. 

3 Reconstruction of the MLDF 



One of the first MLDF parametrizations was introduced by K. Greisen [25], 

-0.75 / ^ \ -2.5 



P,{r) = Mt) ( - 
\r 



1 + ^ 



(9) 



where r is the distance to the shower axis, ro = 320 m and Nfj^{t) is a normal- 
ization constant that depends on the atmospheric depth t. Subsequently other 
groups proposed different functional forms [26,27]. Although all these formula 
describe the MLDF very accurately in the range of short and intermediate dis- 
tances, they are not so good at larger distances and for higher-energy show- 
ers. Recently the KASCADE-Grande Collaboration proposed a new MLDF 
parametrization which is a modification of the Greisen formula [28], 



r 



1 



10 ro. 



(10) 



where A^^, ro, a, Pfj, and 7 are parameters which define the shape and size of 
the MLDF. 
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Fig. 5. Lateral distribution functions of muons at 2.5 m underground for showers 
initiated by protons and iron nuclei of ^ = 30° and primary energies 
and 10^^'^ eV. To generate the air showers we used Aires 2.8.2 with QGSJET-II. 
The solid lines correspond to fits with the KASCADE-Grande MLDF. 

To study the shape of the underground MLDF we performed air shower sim- 
ulations. Again, we used Aires 2.8.2 with QGSJET-II and to propagate the 
muons in the soil we used Eq. (1). Figure 5 shows the simulated MLDFs at 
2.5 m underground corresponding to showers initiated by protons and iron 
nuclei of 30° of zenith angle and for several primary energies. It also shows 
the fits of the profiles with the KASCADE-Grande MLDF where we fixed the 
parameters Tq = 320 m and a = 0.75 and left iV^, and 7 as free fit parame- 
ters. The simulated profiles are fitted very accurately in the range of distances 
considered. 

The range of primary energies for which we simulated MLDFs in order to fit 
them with the KASCADE-Grande formula was 10^^-^ eV - 10^^'^ eV in steps 
of A\og{E/eV) = 0.1. This procedure was performed for both protons and 
iron nuclei at ^ = 30° and ^ = 45°. 

The parameter 7 in Eq. (10) describes the behavior of the MLDF at relatively 
large distances from the core, where only low statistics are available in general. 
On the other hand, (3/^ is well sampled by the proposed 750 m-array. 

Therefore, for all practical purposes, the analysis shows that a better fit of 
the muon counters data to the complete MLDF function is obtained by just 
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fixing 7 = 3, its average value for protons, while leaving and A^^ as free 
parameters. The factor involving 7 only becomes important at distances to 
the shower core larger than ~ 10 x tq = 3.2 km, well beyond the point where 
data is available at the energies of interest. To keep 7 free when dealing with 
real data would only add degrees of freedom to the fitting function without 
the corresponding increase in the available data. 

To develop the reconstruction method it is necessary to simulate an array of 
Cherenkov detectors with associated muon counters. The Cherenkov detectors 
give the trigger information, geometry and energy reconstructions, synchro- 
nization, and telecommunications. For this purpose we interfaced an ad-hoc 
developed muon counter computer code with the program SDSim version v3r0 
[32], which simulates the Pierre Auger surface detectors. 

Since the number of particles produced in a shower is extremely large, ~ 
10^ in a 10^^'^ eV shower, it is practically impossible to follow and store all 
the information of the secondary particles. Therefore, a statistical method 
called thinning, first introduced by M. Hillas [33,34], is used. An unthinning 
method has to be applied to calculate the real number of muons arriving at 
a counter. We employed the method introduced by P. Billoir described in 
references [35,36]. 

The muons with enough energy to reach the detector that arrive in a time 
interval At are uniformly distributed and those that fall into the same segment 
are counted as one, in this way we include in the simulation the pile-up effect 
introduced by the electronics. The present simulation does not include the 
propagation of the electromagnetic particles into the soil. Detailed numerical 
simulations involving Geant4 [37] and Aires show that, beyond 17 radiation 
lengths, the electromagnetic contribution is at most a few percent [38] and 
can therefore be neglected in the present analysis. Nevertheless, these punch 
through simulations will be experimentally verified by a specific detector [39] . 

We used the generated showers to study the shape of the MLDF to simulate 
the response of the Cherenkov detectors and muon counters by randomly dis- 
tributing impact points in the 750 m-array. We used muon counters of 30 m^, 
192 segments and a time resolution of 20 ns. For each muon counter we cal- 
culated the "real" muon distribution (i.e., assuming an infinite segmentation, 
without pile-up cfi^cct), the "measured" muon distribution (taking into account 
the pile-up effect) and the corrected muon distribution (obtained applying Eq. 
(6) in each time bin). 

As an example, figure 6 shows the "measured" time distribution of muons in 
a muon detector at 227 m from the shower core for a simulated event initiated 
by an iron primary of £^ = 10^^ eV and 6 = 30°. Figure 6 also shows the "real" 
time distribution and the corrected one. 



12 



80 



70: 
60: 
50: 
40 ; 

30 ; 

20 



Measured 
I I Corrected 



E = 10^^eV 



9 = 30° 
r = 227m 



10 



III' I I I I 1 1 1 I 1 1 1 I I 



50 



100 



150 
t[ns] 



200 



250 



300 



Fig. 6. Real, measured and corrected time distribution of muons for the closest 
station corresponding to one simulated event oi E = 10^^ eV and 6^ = 30°. We 
considered muon detectors of 30 m^ of area, 192 segments and a time resolution of 
20 ns. 

In this example, the total number of muons corresponding to the real distri- 
bution is ra^^"^ = 172 and the measured one is n^^'^ = 147, 15% less approxi- 
mately. The total number of muons after applying the correction is n'^^^^ = 163, 
about 5% less than the real value. This shows the importance of the correction 
in the determination of the number of muons at each counter, especially at 
those close to the shower core. 

Due to the pile-up effect, counters are not able to measure the number of muons 
very close to the shower axis. Therefore and to avoid systematic uncertainties 
it is convenient to fit the MLDFs including stations with a large number of 
muons in a different way. We will consider as saturated stations those with 
> 72 (i.e. N^'"''' > 90) in at least one time bin. 

The probability of having more than 3 background triggers due to random 
coincidences in a time window At is, P{n > 3) = 1 — exp(— /i)(l + fi + /i^/2), 
where = Iq -^mc At, Iq is the vertical intensity of the background and A^c is 
the muon counters area. The background flux at ~ 2.5 m underground in the 
site where the muon counters are going to be installed is under study. However, 
if we assume a vertical intensity like the corresponding to the Auger tanks (see 
Ref. [40]) we obtain P{n > 3) ~ 8 x 10"'' and for a vertical intensity 10 times 
larger, which can be considered as an upper limit, P{n > 3) ~ 7 x 10""^. 
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Therefore, stations with 0, 1 or 2 muons will be considered as silent stations 
to prevent errors coming from such random coincidences. We will take into 
account all silent stations with distance to the shower core less than 5000 m. 

The number of muons that hit a given detector follows a Poisson distribution. 
Therefore, to fit the MLDF we minimize the likelihood function C = — ln(P) 
with respect to the parameters p= (A^^, of the KASCADE-Grande MLDF 
(see Eq. (10)), where 



ljcxp(-p^(r,;p)) X 

1=1 'Vi ■ 



X 



1 



n exp(-p^(rj;p)) (l + p^{n]p) + -pf^iu^p)'^ j . (11) 

Here the first factor corresponds to saturated stations, the second to "good" 
stations (neither saturated nor silent) and the third to silent ones, is the 
distance of the i-th station to the shower core, Nsat is the number of saturated 
stations, N is the number of "good" stations, Nsu is the number of silent 
stations, is the total number of muons measured corresponding to the 

i-th station and n'^°l^ is the total number of muons for the i-th station after 
applying the correction in every time bin of the measured time distribution. 

The expression for the saturated stations comes from the fact that the Poisson 
distribution of mean value p can be approximated by a Gaussian distribution 
of mean value p and a — ^/p, P{n; p) = exp{—{n — pY /2p) /^/2¥p. Therefore, 
the probability that n be greater than a given no is. 



1 °° 

P{n > Uq) = / o?nexp 

\/27ru J 



{n-pf\ 1 



J27rpj \ 2p 2 

no ^ ' 



1-Erf 



,(12) 



where. 



Erf(a;) = 4= J dtexp{-t^). (13) 

^ 



The expression for the silent stations corresponds to the probability that the 
number of muons be less or equal than two, 

P{n <2) = j2 exp{-p) ^ = exp(-/.)(l + // + ^). (14) 
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Fig. 7. Fit of the MLDF for the same event considered in Fig. 6. The arrows indicates 
the position of the silent stations. We considered muon detectors of 30 of area, 
192 segments and a time resolution of 20 ns. 

Figure 7 shows a fit of the MLDF for the same event of figure 6. To reconstruct 
the arrival direction and the position of the shower core we used the standard 
package, called CDAS [41], used to reconstruct the information from the Auger 
Cherenkov detectors. The arrows correspond to the positions of the silent 
stations. 



In order to stress the importance of including the silent and saturated stations 
in the fit procedure, figure 8 shows the fits corresponding to proton events 
of E = 10^^ eV and 6 = 30°, with (solid line) and without (dashed line) 
taking into account the silent and saturated stations. We observe that the 
fitted MLDF is much more similar, specially close to 600 m from the shower 
axis, to the real one (histogram of the figure) for the case in which the silent 
and saturated stations are included in the fit. The exclusion of the silent and 
saturated stations produces a spurious fiattening of the fitted MLDF, specially 
for those cases in which the closest station is saturated, as clearly shown in 
bottom panel of figure 8. The discrepancies between the fitted MLDF obtained 
without including the silent and saturated stations and the real one increases 
for decreasing primary energies. It introduces important biases in particular 
in the iV^(600) distribution. The inclusion of the saturated and silent stations 
avoid these biases and reduce the error in the determination of iV^(600). 
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Fig. 8. Fit of the MLDFs for a proton shower of E = 10 eV and 9 = 30° injected in 
two different positions of the array. We considered muon detectors of 30 of area, 
192 segments and a time resolution of 20 ns. The arrows indicate the position of the 
saturated (pointing upwards) and silent (pointing downwards) stations. The sohd 
blue line corresponds to the fit with saturated and silent stations and the dashed 
black line corresponds the fit without saturated and silent stations. The event of 
the upper panel does not have saturated stations and the one of the bottom panel 
has both saturated and silent stations. 
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Fig. 9. Relative error of A^^'^^(600) as a function of the logarithm of the energy 
for protons and iron nuclei of 30° (left panel) and 45° (right panel). We considered 
muon detectors of 30 of area and 192 segments. We assumed a time resolution of 
20 ns (solid lines) and a value much greater than the width of the time distribution 
of muons (dashed lines) which corresponds to the most unfavorable case. 

4 Muons at 600 m from the Shower Core 



The number of muons at a given distance from the shower axis has been used 
in the past as a parameter for composition analysis. The AMIGA 750 m infill 
has a detector spacing appropriate to evaluate the number of muons at 600 m 
from the shower core [42]. 

The discrimination power of A^^(600) depends strongly on its reconstruction 
uncertainty. Therefore, to study the uncertainty introduced by the reconstruc- 
tion method we define, 

A^^^'^(600) , , 

ArKea/(600) ' ^ ' 



where A''^'^'^(600) is the reconstructed number of muons from a MLDF fit and 
^Keai^gQQ) is the expected average number of muons from sampling muons in 
a 40 m wide ring (in the shower plane), both at 600 m from the shower axis. 

Figure 9 shows o(e) (relative error of Ai"^^^(600)), obtained from Gaussian fits 
of the distributions of e, as a function of the logarithm of the energy for proton 
and iron primaries and zenith angles of 30° and 45°. We considered muon 
counters of 30 m^, 192 segments and a time resolution of 20 ns (sohd hues). 
As expected, the relative error decreases with energy, due to the fact that 
the number of muons in the showers increases almost linearly with primary 
energy, therefore, the number of triggered muon detectors (with iV^ > 3) also 
increases. 
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From figure 9 we can also see that the relative error for the iron nuclei is smaller 
than the corresponding one for protons because heavier nuclei produce more 
muons. 

To study the robustness of the reconstruction method we assumed a time 
binning much larger than the width of the time distribution of muons, the 
worst case. Figure 9 shows the obtained results (dashed lines), the relative 
error does not increase much, < 3% with respect to the case of 20 ns. 

Note that for the two cases considered the mean value of e, i. e. the bias in 
the determination of N^{600), kept less than ~ 3% in the whole energy range 
considered (from 10^ ''■^ eV to 10^*-^ eV). 

The uncertainty of A^^^'^(600) is less than 20% in the energy range in which the 
750 m-array will be effective, from 10^^'^ eV to 10^^'^ eV. Moreover, it decreases 
with energy, such that for energies greater than 10^^ eV it is smaller than 13%. 
The error in the reconstruction of the primary energy is typically around 
20% [43], and because the number of muons increases almost linearly with 
the energy this uncertainty will dominate in the determination of A^^^"^(600). 
Therefore, muon counters of 30 m^ and 192 segments are acceptable values for 
the design parameters. 



5 Pcirameters Sensitive to the Chemical Composition 

Several parameters obtained from the surface and fluorescence detectors are 
used for the identification of the primary. The difference between parameters 
corresponding to different primaries is due to the fact that showers initiated 
by heavier primaries develop earlier and faster in the atmosphere and also 
have a larger muon content. 

From the distributions of A''^^'^(600) corresponding to a given energy, type of 
primary and zenith angle, we calculated the mean value (iV^"^(600)) and the 
regions of 68% and 95% of probability. Figure 10 shows (Ai'^^'^(600)) and the 
regions of 68% and 95% probabihty (parameter P in the figure) as a function 
of the energy for protons and iron nuclei of 30° of zenith angle. Although the 
distributions for iron and proton show an overlap, the separation is good at 
68% confidence level. 

Besides the muon content of the shower and the depth of the maximum there 
are other parameters that are used for composition analysis [44]. These are the 
structure of the shower front (in particular the rise-time of the signals in the 
surface detectors), the radius of curvature, R, and the slope, (3, of the lateral 
distribution function of the total signal deposited in the water Cherenkov 
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log(E/eV) Iog(E/eV) 

Fig. 10. log(iV^'^'^(600)) as a function of log(S/eV) for protons and iron nuclei of 
6 = 30°. P gives the probability that an event falls in the corresponding shadowed 
region. 

detectors. 

For any given event, we can define a parameter related to the shower front 
structure, using the rise-time of the signals in a selected subset of the triggered 
water Cherenkov detectors, 

1 , , , /400 m\2 ^ ^ 

ti/2 = ^ E(4o - t\o) X j , (16) 



where Nt is the number of stations with signal greater than 10 VEM 0, t\Q 
and are the times at which 10% and 50% of the total signal is collected, 
respectively, and is the distance of the i-th station to the shower axis. Only 
stations at a distance to the shower axis greater than 400 m are included in 
Eq. (16). 

To obtain X^ax, including the effect of the detectors and the reconstruction 
procedure, we assumed that the distribution of the reconstructed Xmax is a 
Gaussian with an energy dependent a as given in Ref. [20]. Therefore, we ob- 
tained the distributions of the reconstructed X^ax from the simulated showers 
by sampling a value from a Gaussian distribution with the mean value given 
by Xmax calculated internally in Aires and a from the interpolation of the 
simulated data of Ref. [20] . 

To compare the discrimination power of the different parameters considered 

^ Vertical Equivalent Muon, signal deposited in a water Cherenkov tank when fully 
traversed by a muon vertically impinging in the center of the tank [40] 
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Fig. 11. Linear fits of r/ as a function of log(£'/eV) for the different parameters 
considered and for = 30° and 9 = 45°. 



{q = N^{600),Xmax, P, R, Sqoo) we define, 



r]{q) 



(17) 



where qA is the parameter for the nucleus A, {qA} is the mean value and a{qA) 
the standard deviation for the distribution of qA- From the definition of 77(g) 
we see that the larger its value, the greater the discrimination power of the 
parameter q. 

Presumably, the interpolated signal of the water Cherenkov detectors at 600 
m from the shower axis 6*600 will be used to obtain the primary energy for 
the 750 m-array. Therefore, to study its dependence on the primary mass we 
included it in the set of parameters. Figure 11 shows the linear fits of 17 as a 
function of the logarithm of the energy for 9 = 30° and 6 = 45° and for the 
different parameters considered. 

From figure 11 we see that the parameter which better separates protons from 
iron nuclei appears to be Nl^'^^{600), followed by Xmax- We can also see that, 
except for Xmax, V increases with the primary energy. This happens because, 
as the energy increases, the number of triggered stations also increases, which 
reduces the reconstruction errors. The number of particles in the detectors 
also increases with the energy and therefore reduces fluctuations. 

Although the reconstruction error of X^ax decreases with primary energy, rj 
slowly decreases. This happens because the mean values of X^ax corresponding 
to protons and iron nuclei get closer as the energy increases for the hadronic 
model considered (however, this is not the case for the hadronic model EPOS, 
see Ref. [45]). 
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Fig. 12. ry as a function of log(£'/eV) for the parameters Xjnax and iV|^'^^(600) for 
6 = 30° and 9 = 45°, assuming a power law spectrum of spectral index 7 = 2.7 and 
a gaussian uncertainty of 20% of relative error and without including the effects of 
the energy uncertainty (see appendix A for details). The differences between the 
curves corresponding to X^^ax with and without including the energy uncertainty 
are so small that cannot be seen in both plots. 

If the energy uncertainty is negUgible, iV^'^'^(600) is the best parameter for 
mass discrimination analysis. However, as opposed to the other parameters 
which depend on the energy logarithmically, the number of muons is almost 
proportional to the primary energy. Therefore, the energy uncertainty will 
affect more the discrimination power of Nl^^%600) than that of the other 
parameters. Figure 12 shows 77 as a function of the energy assuming a 20% 
energy uncertainty for the parameters iV^'^'^(600) and Xj^ax (see appendix A 
for the details of the calculation). We see that the discrimination power of 
jyAec^-ggg^ decreases in such a way that it is of the order of the corresponding 
to Xmax- From figure 12 we also see that the discrimination power of X„iax 
remains approximately the same when the energy uncertainty is included, 
which is due to its logarithmic dependence on the primary energy. 



6 Conclusions 



In this paper we studied the effect of muon counter segmentation and we 
found an analytical expression to correct for the undcrcounting due to muon 
pile-up. We presented a detailed method to reconstruct the muon lateral dis- 
tribution function for the 750 m infill array. We also studied the potential 
of the A^f^'^(600) parameter as a discriminator of the identity of the primary 
cosmic ray particle. 

We showed that for 192 segments and 90 incident muons in a 20 ns time bin 
(maximum number of muons at 200 m from the shower axis for E = 10^^ 
eV, which would correspond to saturation for the adopted configuration) the 
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total error, Poissonian plus pile-up correction, does not exceed the Poissonian 
error by more than 14%. This correction is enough to attain, with surface 
information alone, a composition discrimination power comparable, or even 
better, than that of Xj^ax- Clearly, this result depends quantitatively on the 
assumed hadronic interaction model and, consequently, the safest conclusion 
is probably that a two-dimensional analysis, including both parameters, might 
well configure an optimal composition analysis method which will be discussed 
in an accompanying paper [46]. 

The error in the determination of A^^^'^(600) parameter is < 20% at lO^'^-^ eV 
(energy threshold of the 750 m array) and decreases as the energy increases, 
being less than 12% for E > 10^^ eV. So, at this energy or greater, this 
uncertainty is small compared with the expected primary energy uncertainty, 
~ 20%, which affects directly the determination of A^f'^'^(600) since shower 
muon content is nearly proportional to the primary energy. Therefore, 30m^ 
individual muon detectors partitioned in 192 segments, can be considered as 
an optimal compromise between cost and scientific output. 

We also studied the effect of the energy uncertainty on the discrimination 
power of A^f^'^(600). The inclusion of this effect somewhat degrades the dis- 
crimination power of A^^^'^(600), however, it generally remains as good a pa- 
rameter as Xmax- 
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A Effect of the Energy Uncertainty 

We estimate the effect of the energy uncertainty on the discrimination power 
of the parameters A^^^'^(600) and X^ax, assuming a power law spectrum of 
spectral index 7 = 2.7, J{E) oc E~'^, and a gaussian uncertainty in the energy 
determination of a = sE with e the relative error. Let g be a parameter for 
which we know its distribution function parametrized by the primary energy, 
f{q;E). Therefore, the conditional probability of q given the reconstructed 
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energy Erec is, 

E2 



P{q \Erec) = C{Erec) J dE f{q; E) exp E-^-\ (A.l) 

where Ei = 10^^ eV and E2 = 10^° eV are the lower and upper hmits of the 
part of the spectrum considered, respectively, and, 

C-\Erec) = JdE exp E-'-\ (A.2) 

El ^ ^ / 



From Eq. (A.l) we can calculate the mean value and the variance of the 
parameter g as a function of the reconstructed energy, 



{q){E,e) = C{Erec) JdE {qom exp P^^^^/^^ ) (^-3) 

El 
E2 

Var[q]{E,e)^C{Erec) J dE iVar[qo]iE) + {qo)\E)) x 

El 

exp E-^-' - mE,e), (A.4) 

where {qo){E) and Var[qo]{E) are the mean value and the variance of param- 
eter q without taking into account the energy uncertainty, calculated using 
f{q;E). Therefore, to calculate the parameter 1] defined in Eq. (17) includ- 
ing the energy uncertainty, we need the functions {qo){E) and Var[qo]{E) for 
each parameter considered. For that purpose, we fitted the mean value and 
the standard deviation of A^^^'^(600) with a function of the form g{E) — a E^ 
and for X^ax with h{E) = a + &log(£'). Using the fits obtained for 9 = 30° 
and 9 = 45° and Eqs. (A. 3), (A.4) and (17) we obtain r]{Erec,^) which can be 
seen in Figure 12 for e = 0.2. 
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